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So far, Brillouin scattering and cavity optomechanics were mostly disconnected branches of re¬ 
search - although both deal with photon-phonon coupling. This begs for the development of a 
broader theory that contains both fields. Here, we derive the dynamics of optomechanical cavi¬ 
ties from that of Brillouin-active waveguides. This explicit transition elucidates the link between 
phenomena such as Brillouin amplification and electromagnetically induced transparency. It proves 
that effects familiar from cavity optomechanics all have traveling-wave partners, but not vice versa. 
We reveal a close connection between two parameters of central importance in these fields: the Bril¬ 
louin gain coefficient and the zero-point optomechanical coupling rate. This enables comparisons 
between systems as diverse as ultracold atom clouds, plasmonic Raman cavities and nanoscale sili¬ 
con waveguides. In addition, back-of-the-envelope calculations show that unobserved effects, such as 
photon-assisted amplification of traveling phonons, are now accessible in existing systems. Finally, 
we formulate both circuit- and cavity-oriented optomechanics in terms of vacuum coupling rates, 
cooperativities and gain coefficients, thus reflecting the similarities in the underlying physics. 


I. INTRODUCTION 

Brillouin scattering [1] and cavity optomechanics [2] 
have been intensively studied in recent years. Both con¬ 
cern the interaction between light and sound, but they 
were part of separate traditions. Already in the early 
1920s, diffraction of light by sound was studied by Leon 
Brillouin. Therefore, such inelastic scattering is called 
Brillouin scattering [3, 4]. The effect is known as stim¬ 
ulated Brillouin scattering (SBS) [5-7] when a strong 
intensity-modulated light field generates the sound, of¬ 
ten with classical applications such as spectral purifica¬ 
tion [8] and microwave signal processing [9] in mind. In 
contrast, cavity optomechanics arose from Braginsky’s 
efforts to understand the limits of gravitational wave de¬ 
tectors in the 1970s - and greatly expanded since the 
demonstration of phonon lasing in microtoroids [10]. By 
and large, it aims to control both optical and mechanical 
quantum states [11-13]. 

Historically, a number of important differences hin¬ 
dered their merger. For instance, SBS generally dealt 
with high-group-velocity and cavity optomechanics with 
low-group-velocity acoustic phonons. In addition, bulk 
electrostrictive forces usually dominated phonon gener¬ 
ation in SBS - while radiation pressure at the bound¬ 
aries took this role in cavity optomechanics. Further, 
cavity optomechanics typically studied resonators with 
much lower phonon than photon dissipation - whereas 
Brillouin lasers [8, 14, 15] operate in the reversed regime 
[16]. Finally, SBS is often studied not in cavities but 
in optically broadband waveguides [1]. Thus, particular 
physical systems used to be firmly placed in either one 
or the other research paradigm. 

Lately, the idea that these are mostly superficial classi¬ 
fications has been gaining traction. Indeed, in both cases 
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FIG. 1 . From circuit to cavity optomechanics. We ex¬ 
plicitly derive the physics of optomechanical cavities (right) 
from that of Brillouin-active waveguides (left). Therefore, 
both traveling-wave and cavity-based optomechanics can be 
cast in terms of vacuum coupling rates (go and go), coopera¬ 
tivities ( C and C) and gain coefficients (Q and Q). 

light generates motion and the motion phase-modulates 
light. Next, this spatiotemporal phase-modulation cre¬ 
ates motional sidebands - which interfere with those 
initially present. The research fields share this essen¬ 
tial feedback loop. Some connections have already been 
made. For instance, electrostrictive forces were exploited 
for sideband cooling [17, 18] and induced transparency 
[19, 20] while radiation pressure contributed to SBS in 
dual-web fibers [21] and silicon waveguides [22-25]. 

In this work, we derive the dynamics of optomechan¬ 
ical cavities from that of Brillouin-active waveguides 
(fig.l). The transition holds for both co- and counter- 
propagating pump and Stokes waves, i.e. for “for¬ 
ward” and “backward” scattering, and for opto-acoustic 
coupling between two different or two identical optical 
modal fields, i.e. for “inter-” [18, 26-29] or “intra-” 
[22, 23, 28, 30-35] modal scattering (fig.2). Hence, all 
flavours of photon-phonon interaction are treated on the 
same footing. Moreover, this spatially averaged cav¬ 
ity dynamics is found to be equivalent to the standard 
Hamiltonian of cavity optomechanics [2] - even in the 
case of low-finesse phonons. It turns out that this cav¬ 
ity dynamics can be mapped - by swapping space and 
time (z t) - on the steady-state spatial evolution of 
the opto-acoustic fields in the waveguide. The treatment 
describes Raman effects as well. 
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FIG. 2. Phase-matching diagrams. The optical disper¬ 
sion relation ui(k) shows phonon-mediated coupling between 
co- or counter-propagating photons and between two identical 
(intra-) or two different (inter-modal) optical modes. There¬ 
fore, there are generally four types of optomechanics, of which 
three are indicated in this diagram. The fourth is counter¬ 
coupling between two different optical modes. 
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FIG. 3. Symmetry of circuit and cavity optomechan¬ 
ics. Each temporal optomechanical effect has a spatial sym¬ 
metry partner. Thus, the description of these effects can be 
cast in terms of conceptually similar figures of merit. The 
scheme assumes a red-detuned optical probe; “gain” and 
“laser” should respectively be replaced by “loss” and “cool¬ 
ing” for a blue-detuned optical probe. The meaning of the 
figures of merit is discussed in the main text. 


This implies that the plethora of optomechanical ef¬ 
fects, such as stimulated Brillouin scattering [3, 4, 
36], slow light [37-39], optomechanically induced trans¬ 
parency [39, 40], ground-state cooling [11, 41] etc., are 
different aspects of the same feedback loop. The rig¬ 
orous transition decisively indicates that both fields are 
a subset of a larger theory of photon-phonon interac¬ 
tion, which may be built on a single Hamiltonian [42]. 
This is not to say that they are identical: a Brillouin- 
active waveguide supports complex spatiotemporal phe¬ 
nomena [43-45] and noise dynamics [46, 47] not present 
in a high-finesse optomechanical cavity. Nevertheless, 
in the resulting picture (fig.3), both traveling-wave and 
cavity-based photon-phonon interaction can be classified 
according to (1) the damping hierarchy of the photons 
and phonons and (2) the strength of the photon-phonon 
coupling with respect to the largest dissipation channel. 
For weak coupling, the long-lived particle species - ei¬ 
ther photons or phonons - triggers the photon-phonon 
conversion. The short-lived particle species cannot truly 
build up and is thus slaved to its long-lived partner; it is 
merely created in short segments (of space or time) and 
immediately decays afterwards. 

All Brillouin-active waveguides so far exhibited far 
stronger phononic than photonic propagation losses; 
in addition, the coupling was always weak relative 


to this phononic damping. Hence, there are two to 
date unexplored regimes of guided-wave optomechanics: 
(1) photon-assisted amplification of traveling phonons 
and (2) strong coupling between traveling photons and 
phonons (fig.3). The strong coupling regime produces 
either traveling entangled photon-phonon pairs or state 
swapping between light and sound along the waveguide, 
depending on the details (e.g. probe detuning) of the 
experiment. Although currently unobserved, both ef¬ 
fects may be an asset in future quantum phononic net¬ 
works [13, 48-51]. For instance, in the strong coupling 
regime the flying phonon - entangled to its photonic 
partner - could be detected piezo-electrically or optically 
and thereby enable Bell tests [52-54] between two differ¬ 
ent particle species. Our back-of-the-envelope estimates 
show that these regimes can be achieved in existing sys¬ 
tems, such as dual-web fibers and silicon nanowaveguides. 

The transition (fig. 1) assumes that the photonic and 
phononic modes of the waveguide are not disturbed too 
strongly by looping it into a cavity. This is justified 
in many cases since cavity designs aim to minimize the 
losses (e.g. due to bending) induced by any modal per¬ 
turbations. Within this approximation, it permits trans¬ 
lations between circuit- and cavity-oriented figures of 
merit. For instance, we identify a connection between 
the Brillouin or Raman gain coefficient Q and the zero- 
point coupling rate go- The former (Q) quantifies the 
pump power and waveguide length required to amplify a 
Stokes seed appreciably [3, 4], The latter (go) captures 
the interaction strength between a single photon and a 
single phonon in an optomechanical cavity [2]. The tran¬ 
sition proves that these figures of merit are inextricably 
linked by 
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with v p and v s the group velocities of the pump and 
Stokes waves, hu > p the pump photon energy, the 
mechanical resonance frequency, L the cavity roundtrip 
length and Q m the waveguide’s mechanical quality fac¬ 
tor. This link is independent of the type of driving op¬ 
tical force and of the relative photon and phonon damp¬ 
ing. Similarly, we derive connections between each of the 
circuit- and cavity-oriented figures of merit: between the 
vacuum coupling rates (go and go, see (21)), the cooper- 
ativities (C and C, see (34)) and the gain coefficients (Q 
and Q , see (36)). 

Notably, this treatment goes beyond cavity optome¬ 
chanical systems that have a clear circuit equivalent 
(as in fig. 1). Indeed, the standard cavity Hamiltonian 
H = huj c (x)a^a + Ml m b'b [2] also captures the tempo¬ 
ral dynamics of cavity optomechanics based on Bose- 
Einstein condensates [55, 56] or plasmonic Raman cavi¬ 
ties [57]. The physics of all these diverse systems can be 
understood in the scheme of fig.3. On top of the similar 
dynamics, this means that the photon-phonon interac¬ 
tion efficiency of a larger class of systems can now be 
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FIG. 4. Example of phonon dispersion relation, a, The 

frequency f1(g) of transversally trapped acoustic phonons has 
a Raman-like cut-off fi c for low q and approaches the bulk 
relation for large q. b, Thus, the phonon group velocity v m 
vanishes for low q and becomes the bulk speed for large q. 

compared in a single framework. For instance, the gain 
coefficient of a silicon nanowire can be converted to the 
vacuum coupling rate of a hypothetical cavity (through 
( 1 )); which can next be compared to that of any other 
cavity optomechanical system. In reverse, the link en¬ 
ables the conversion of a vacuum coupling rate of an 
actual cavity optomechanical system into a hypotheti¬ 
cal guided-wave coupling rate (through ( 21 )); which can 
next be compared to that of any other waveguide. We 
give examples of such conversions, which can be tested 
empirically in many cases, in section V. 

The paper is organized as follows: in section II we de¬ 
scribe a minimal model of traveling-wave optomechanics 
and frame it in terms of a guided-wave vacuum coupling 
rate go and cooperativity C. Next, we make the mean- 
field transition to a cavity in section III. At that point, 
we also discuss the limitations of the analysis. The re¬ 
sulting dynamical effects are treated in section IV. The 
prospects for observing new effects are considered in sec¬ 
tion V and we conclude in section VI. 


II. CIRCUIT OPTOMECHANICS 

In particular, we study the interaction between a pump 
field with envelope a p (z , t ) and a red-detuned Stokes field 
envelope a s (z, t ) mediated by an acoustic field with enve¬ 
lope b(z,t). The envelopes contain only the slowly vary¬ 
ing part of the photonic-phononic fields; rapidly oscil¬ 
lating factors e *( fez -^) were removed in each case. The 
guided optical modes correspond to the points (w p ,fc p ) 
and (w s ,fc s ) in the optical dispersion relation (fig.2). By 
energy and wave-momentum [58] conservation, the ex¬ 
cited phonon has an angular frequency ft = ui p — cj s and 
wavevector q = k p =p k s . The nature of the optical modes 
(co/counter and fast/slow) and the acoustic dispersion 
relation determine the wavevector q and group velocity 
v m of the excited phonons (fig.2&4). 

Traveling-wave photon-phonon coupling is governed by 


the following dynamical evolution [42, 47, 59] 

v~ 1 d t a p + d z a p = -ig 0 a s b - Xp 1 ^ 
v~ 1 d t a s ± d z a s = -igotf a p - x” 1 ^ ( 2 ) 

v^dtb + d z b = -ig 0 ala p - 

Its derivation proceeds from Maxwell’s and the elastic¬ 
ity equations on the assumption that the envelopes vary 
slowly in space and time [42, 47, 59]. This starting point 
and the following treatment holds quantum mechanically 
if one takes care to treat the envelopes in ( 2 ) as operators 
[42, 47] obeying the equal-time commutator 

[ffi 7 (z,f),aj y ,(z , ,t)] = ^/VyVy5jy5(z - z') (3) 

with “ 7 ” an index running over the pump “p”, Stokes 
“s” and mechanical wave “m”, u 7 the group velocities, 
<5 7 y the Kronecker delta, 8{z) the Dirac delta distribu¬ 
tion and a m = b for notational convenience. We flux- 
normalized the field operators a 7 such that <h p = aj,a p , 
<f> s = a|a s and <F m = b^b correspond to the number 
of pump photons, Stokes photons and phonons passing 
through a cross-section of the waveguide per second. We 
will treat highly occupied (i.e. large mean flux (4> 7 )) 
modes as classical in the remainder of the paper, as is 
standard [2, 60-62]. Further, we denote go the traveling- 
wave vacuum coupling rate (to be discussed further on), 
X 7 1 = ^2 ~ *A 7 the susceptibilities, a 7 the propagation 
losses and A 7 the wavevector offsets between externally 
applied fields and the intrinsic waveguide modes. 

In some systems, e.g. for the Raman-like low-group- 
velocity phonons (fig.4) associated with forward intra- 
modal scattering [22-24, 28, 30], the phonon wavelength 
can be substantially larger than its decay length a " 1 - 
so its slowly-varying amplitude treatment breaks down. 
Then the acoustic excitation is better treated as a lo¬ 
calized series of mechanical oscillators [23, 30, 59], essen¬ 
tially dealing with each cross-sectional slice of the waveg¬ 
uide as an artificial Raman-active molecule. The above 
dynamics ( 2 ), however, contains these systems as well by 
letting the phonon decay length a " 1 vanish. Further, the 
sign (±) in the Stokes equation indicates the difference 
between forward (+) and backward (—) photon-phonon 
coupling. Cascaded scattering [30, 63] and noise [46, 47] 
can and should be added to this model in some instances. 
In fact, (2) can be regarded as the unique, minimal model 
for guided-wave Brillouin scattering [42, 47, 59]. We dis¬ 
cuss potential extensions in section V; in the following, 
we only need the minimal model ( 2 ), future extended 
versions can be dealt with similarly. 

The Manley-Rowe relations [3, 64] guarantee that a 
single, unique figure of merit go captures all conservative 
optical forces and scattering. Indeed, in the lossless case 
(a 7 = 0 ), the rate of pump photon destruction must 
equal the rate of Stokes photon and phonon creation: 

-<9 2 $p = ±<9 2 <F s = d z $ m = -% (ialtfa p + h.c.) (4) 
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Similar to go in a cavity [2], go quantifies the inter¬ 
action strength between a single photon and a single 
phonon, but in this case flying along a waveguide in¬ 
stead of trapped in a cavity. We take go real and positive 
without loss of generality. Briefly specializing to forward 
intra-modal scattering, the mean-field transition of sec¬ 
tion III will show that (see Appendices) 


with 


go = -xzpf 


dk p 

dx 


( 5 ) 


a?ZPF = £zpf 



h 




( 6 ) 


the guided-wave zero-point motion and m e ff the effec¬ 
tive mass of the mechanical mode per unit length. In¬ 
deed, a short waveguide section of length 8L contains 
(n m ) = — (<E> m ) phonons with (<k m ) the mean phonon 
flux. As particle fluxes - instead of numbers - are funda¬ 
mental in the traveling-wave Manley-Rowe relations (4), 
the zero-point motion is rescaled by precisely this factor 
(dL/vm ) 1 / 2 relative to the actual zero-point motion .tzpf 
[2] of the 8L -section 


XZPF - /2 (7) 

Therefore, the traveling-wave vacuum coupling rate go is 
determined by the wavevector shift induced by mechan¬ 
ical motion at fixed frequency, while the cavity vacuum 
coupling rate go is determined by the frequency shift in¬ 
duced by mechanical motion at fixed wavevector [2]. No¬ 
tably, the interpretation of go as the coupling strength be¬ 
tween a single traveling photon and phonon holds also for 
inter-modal and backward scattering (see Appendices). 

In steady-state (<9 t —» 0) and for a constant, strong 
pump ($ p (,z) = $ p ( 0 )), the evolution ( 2 ) reduces to 

&z a s = + X s (8) 

d z b= -ig 0 a\a p - x^b 

The phonon decay length a " 1 is generally largest for 
backward scattering. Even then, it typically does not 
exceed a " 1 ~ 100 /im [3, 65]. Therefore, the photon de¬ 
cay length massively exceeds the phonon decay length 
in Brillouin-active waveguides to date (a s <C a m ). A 
full solution of ( 8 ) exists but yields little intuitive in¬ 
sight (see Appendices). Therefore, we initially focus on 
two subcases: the conventional (a s -C a m ) and the re¬ 
versed case (cc m <C a a ), both in the weak coupling regime 
G?o-\/^p cc s + a m ). These examples illustrate how one 
can formulate guided-wave optomechanics, including the 
classical stimulated Brillouin regime, in terms of the vac¬ 
uum coupling rate go and cooperativity C. 

First, strongly damped phonons (a s <C a m ) act as a 
localized slave wave ( d z b —> 0 ) given by b = —iXmEloatcip- 
On resonance (A 7 = 0), we thus have 

d z a s = +(1 -C) ya s 


with 


£ _ 4gg$p = 4g 2 


( 10 ) 


the guided-wave cooperativity and g = go\f^p the 
pump-enhanced spatial coupling rate. Therefore, C = 1 
is the threshold for net phonon-assisted gain on flying 
photons. Since P p = boj p ^ p is the pump power, we ob¬ 
tain C = ^ and 

O' s 


g = 4 ' g2 ° 

buj p a m 


( 11 ) 


the well-known Brillouin gain coefficient [3, 4], here 
framed in terms of a spatial coupling rate go and cooper¬ 
ativity C. It characterizes the spatial exponential build¬ 
up of a Stokes seed in case of highly damped phonons 
(a s < a m ). Since ($ rn ) m gN?($ s ) <C ($ s ), there are 
on average far fewer phonons than photons flying along 
the waveguide in this case. The system enters the strong 
coupling regime as soon as C ~ ^ (see section IV). 

Second, when the phononic damping is lowest (a m -C 
a s ), we similarly get a slaved Stokes wave ( d z a s —» 0) 
given by a s = —ixsgob^a p resulting in (A 7 = 0) 

d z b=-(l-C) a ^b ( 12 ) 


such that C = 1 also yields the threshold for net 
photon-assisted gain on flying phonons. Since (<E> S ) = 
^ S -C(4> m ) <C (‘km)) there are far fewer photons than 
phonons flying along the waveguide in this case. The sys¬ 
tem enters the strong coupling regime as soon as C ~ 

By replacing the undepleted pump with an undepleted, 
strong Stokes mode (g = goV&s), it follows that an anti- 
Stokes seed sees larger loss by a factor (1+C) convention¬ 
ally (cc s <C a m ) and that a guided-wave phonon channel 
can be cooled by a factor (1+C) when it has the low¬ 
est propagation loss (a m <C a s ). An undepleted, strong 
phononic beam (g = goV^m) yields similar coupling be¬ 
tween the pump and Stokes wave. 

The coupling rate g and the cooperativity C respect the 
symmetry between flying photons and phonons, whereas 
the gain coefficient Q ( 11 ) is most relevant in case of 
stronger phonon damping. Therefore, we regard g and 
C as more natural and fundamental figures of merit. It 
is straightforward to extend the above discussion for ab¬ 
sorptive decay of the pump flux, i.e. 4> p (z) = 4> p (0)e _ap;2 
and non-zero wavevector detunings A 7 ^ 0. 

So far, we discussed two subcases of guided-wave Bril¬ 
louin scattering. We treat the strong coupling regime in 
section IV and the full solution in the Appendices. Next, 
we move on to cavity optomechanics via the mean-field 
transition. 


(9) 
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III. BRIDGE TO CAVITY OPTOMECHANICS 

In this section, we transition to an optical cavity - 
made from a Brillouin-active waveguide - of roundtrip 
length L (fig.l). To do so, we introduce the mean-field 
envelope operators 

a(t) = — f a(z,t)dz (13) 

L Jo 

for both the optical (a p / s (f)) and acoustic (5(f)) fields. 
Such mean-field models have found early use in the treat¬ 
ment of fluorescence [ 66 ] and recently also in the context 
of frequency combs [67]. During roundtrip propagation, 
each field obeys dynamics of the form (see ( 2 )) 

v~ 1 d t a +d z a = C - x~ L a (14) 

with £ the nonlinear interaction term. To describe the 
cavity feedback (fig.l), we add the boundary condition 

a( 0 ,f) = Vl — ot! \J\ — ge lv a(L,t) + \/gs(t) (15) 

with a 1 the additional loss fraction along a roundtrip (on 
top of a , such as bending losses), g the fraction of pho¬ 
tons or phonons coupled to an in- or output channel, <p 
the roundtrip phase shift and s(f) the flux-normalized 
envelope of injected photons or phonons. By Taylor- 
expansion of (15), we get 

a(L, t) - a{ 0, t ) « - iip'j a(t ) - s(t ) (16) 

with higher-order terms negligible and a(L , t) ~ a(t ) in 
the high-finesse limit. Low-finesse situations, particularly 
relevant for phonons, are treated further on (see (24)). 
We operate close to the cavity resonance, such that <p <C 
27t. Next, we let (13) operate on (14) and use <9*a = a(t): 

u - 1 a(t) + L~ 1 {a(L, t) - a(0, t)} = C(f) - x _ 1 a(f) (17) 

We insert (16) in (17) and find 

U = v(-x~ 1 a+ (18) 

with x~ L = f — *A the cavity’s photonic or phononic 
response function, k = Ki + k c the total decay rate, K; = 
a ~f aL the intrinsic decay rate, k c = f the coupling rate, 
A = L the detuning and T = — the roundtrip time. 

Next, we multiply (18) by \/T and switch from flux- 
to number-normalized fields (a H> vTa): 

a = vVT ( — y _1 a + s/ftf s (19) 

From here on, n = a^a represents the number of quanta 
in the cavity, while still corresponds to the injected 
photon or phonon flux. The transition from (14) to (19) 
still holds when we replace z H> — z because condition 




FIG. 5. Cavity description. The photonic and phononic 
density of states The mean-field equations (20) de¬ 

scribe coupling between one phononic and either one (a) or 
two (b) photonic resonances. The latter case (b) is most 
power-efficient, although hard to achieve in practice [8]. 


(16) also reverses. Therefore, potential dynamical dif¬ 
ferences between forward and backward scattering dis¬ 
appear in a high-finesse traveling-wave cavity - at least 
in the minimal model (2) of guided-wave optomechanics. 
Comparing (2) to (14), we see that £ tx fg with / and g 
equal to a p / s or b. In the mean-field approximation, we 
assume these envelopes vary little over a roundtrip such 
that fg = fg holds (see Appendices). Finally, we ap¬ 
ply the mean-field (14)-to-(19) transition to (2). Hence, 
an optomechanical cavity - constructed from a Brillouin- 
active waveguide - is governed by 

G-p — igo^sb ,\ p T ^Sp 

«s = -igob ] a p - x" 1 ^ + (20) 

b — igoala p Xm ^ T \Jo~ :cn] .s m 


with 


v p v s v m 


go = go 


( 21 ) 


the well-known temporal zero-point coupling rate [2]. In¬ 
deed, equations (20) are equivalent (see Appendices) to 
the Heisenberg equations of motion resulting from the 
well-known Hamiltonian f-L = hui c {x)a^a + hfl m mb [2]. 
Remarkably, the equivalence holds even for inter-modal 
and backward scattering. The connection (21) between 
the traveling-wave and the cavity-based vacuum coupling 
rates go and go is at the heart of this work: other links 
such as (1) are based on this result. Further, the mean- 
field transition transforms the guided-wave commutator 
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(3) into 

V -j^-Srf [ \ dzdz'S(z-z') 

L Jo Jo 

vZVy 

-j^- S rr> ( 22 ) 

and through rescaling a 7 by y/T\ into 

j^U 7 , j (23) 

thus correctly retrieving the standard harmonic oscillator 
commutators [2]. 

To derive (20), we made the same mean-field transition 
for photons and phonons. In particular, this supposes a 
large phonon finesse 1. Often there is only 

intrinsic phonon loss such that K m = v m a m and thus this 
requires S> 1. In many systems, the phonon decay 
length a" 1 is much shorter than the roundtrip length 
L. Then this phonon high-finesse limit does not hold. 
However, we can completely neglect phonon propagation 
(i d z b —» 0 in (2)) if a m is sufficiently large. The phonons’ 
envelope operator b then obeys 

v7n 1 dtb = -ig 0 ala p - x^b 

Applying (13), multiplying by y/T m and switching from 
flux- to number-normalized envelopes results in 

b = -ig 0 ala p - v m x^b (24) 

where we used (21). Hence, this localized low-phonon- 
finesse approach yields the same result as the previous 
high-finesse limit with s m = 0 (see (20)). Therefore, 
even low-finesse phonons produce the same dynamics as 
is commonly studied in cavity optomechanics [2]. 

Notably, the standard treatment of cavity optomechan¬ 
ics [2] does not consider an explicit space variable: the 
Hamiltonian H performs an implicit spatial average by 
describing the entire object as single mechanical oscilla¬ 
tor, in contrast to the explicit spatial average (13) per¬ 
formed in this work. However, even the implicit average 
in 7i requires low-loss acoustic excitations to set up a 
global mechanical mode self-consistently, precisely as in 
the high-finesse approximation leading to (20). In the lo¬ 
calized, low-finesse phonon approach that generates (24), 
the spatial averaging can still be performed and yields the 
same classical dynamics - but its meaning changes. Now 
{Tm < 1) the acoustic wave is too lossy to set up a global 
mechanical mode for the entire cavity. Instead, the cav¬ 
ity consists of an ensemble of independent Raman-like 
mechanical oscillators. It is no longer possible to address 
phonons circulating in the cavity. 

Finally, we combine (21) and (11). Using a m v m = 
Um/Qm, we obtain result (1) immediately. Note that 
Q in is defined here as the waveguide’s intrinsic phonon 
quality factor, which could be different from the cavity’s 
phonon quality factor if there were e.g. non-negligible 



phonon coupling or bending losses. In case of doubt, it 
is safe to alternatively write (1) as 


VpV s V m 


(Hw p ) a D 


0 = 9o 


Both Q and go are well-established in the study of photon- 
phonon interaction, but they operate on different levels. 
The Planck constant h enters (1) because Q is classical 
while go is inherently quantum mechanical. In addition, 
Q quantifies the combined action of forces and scatter¬ 
ing and contains the phonon loss - while go does not. 
Further, larger L yields a smaller go while Q is length- 
independent. Therefore, g p oc £ q . This mean-field 

derivation is but one way to prove the Q go conver¬ 
sion, other approaches yield the same result (see Appen¬ 
dices). This proof captures all reversible photon-phonon 
coupling mechanisms. 


IV. SYMMETRY BETWEEN CIRCUIT AND 
CAVITY OPTOMECHANICAL EFFECTS 

In this section, we describe both guided-wave and 
cavity-based regimes of photon-phonon coupling. To be¬ 
gin with, we recover and briefly review the known cavity- 
based regimes of photon lasing, phonon lasing and strong 
coupling. Next, we map these regimes on the guided- 
wave spatial evolution of the opto-acoustic fields. The 
mapping unveils two unobserved regimes of guided-wave 
Brillouin scattering. We pay particular attention to the 
strong coupling regime (g a s + a m ). 

Here, we assume zero photon and phonon input flux 
and an undepleted pump. Then (20) reduces to 

a s = —igotfa p - x s _1 a s (26) 

b = -ig 0 ala p - x^b 

These equations treat the photons and phonons identi¬ 
cally. Therefore, every photonic phenomenon must have 
a phononic counterpart and vice versa. Even more, the 
temporal cavity dynamics (26) can be mapped ( 1 1 —> z) on 
the spatial steady-state waveguide evolution (8). Each ef¬ 
fect known from cavities therefore has a waveguide coun¬ 
terpart (but not vice versa as we will see). This also 
implies that the spatial figures of merit have a temporal 
symmetry partner and vice versa; we prepared for this at 
the end of section II by defining a guided-wave vacuum 
coupling rate go and cooperativity C. To clearly expose 
these symmetries, we solve (26); keeping in mind that 
the very same discussion holds spatially for (8). First, 
we decouple (26) and get 

(s +xr )(s + x ") 5(,) = !,25w <27) 

Here, we introduced the pump-enhanced coupling rate 
9 = fi'Ox/^p- Next, we insert the ansatz b oc e 7t in (27) 
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and find two roots 7± 

7± = \ { - (x s “* + Xm 1 ) ± \J (Xs“* - Xm)' +4ff 2 | 

(28) 

In general, these roots strongly mix the photon and 
phonon response: the photon-phonon pair forms a po- 
lariton [12, 56, 61, 68 , 69]. The guided-wave analog of 

(27) is 

ipz ± Xs"*) ipz + Xm 1 ) b{z) = ±9 2 Kz) (29) 

and it can be treated identically. The full spatial and 
temporal dynamics is governed by the general solution 

(28) (see Appendices). However, it is more instructive to 
consider the limiting cases of weak and strong photon- 
phonon interaction relative to the system’s damping. 

First, if the photon-phonon interaction is sufficiently 
weak, i.e. g -C |/c s — /c m |, the two roots in (28) discon¬ 
nect. Usually, the photon and phonon decay rates differ 
significantly. Then there are two scenario’s depending 
on the relative photonic and phononic decay rates. Es¬ 
sentially, the dynamics of the short-lived particle can be 
adiabatically eliminated, although it may still strongly 
modify the response of its long-lived partner. 

In particular, when the phonons decay slowly (/t m -C 
k s ), the photonic response is barely modified: o s —> 0 

and therefore a s = —iXsgob a P to a good approximation. 
However, the phononic response can then dramatically 
change to y ” 1 — g 2 xt- Hence, we recover the spring effect 
(<5f2 m = g 2 ^sxt) and phonon lasing (<5/t m = —2g 2 5Ry*) [ 2 ]. 
At the photon resonance (A s = 0), the phonon linewidth 

equals K m + Sn ni = (1 — C) K m with C = the coop- 
erativity. Therefore, the threshold for meclranical lasing 
is C = 1. This instability was first contemplated by Bra¬ 
ginsky [70] and received further study in systems ranging 
from gram-scale mirrors [71] to optomechanical crystals 
[72, 73]. Since (n s ) = ^“U(n m ) <C (n m ), there are far 
fewer Stokes photons than phonons in the cavity in this 
situation. The system enters the strong coupling regime 
as soon as C ~ ■^ s ~. 

Similarly, when the photons decay slowly (k s <C K m ), 

the phononic response is barely modified: b —> 0 and 
therefore b = — iy m 5 oh|a p to a good approximation. 
However, the photonic response can then dramatically 
change to y^ 1 — 9 2 X m- Hence, we recover the cav¬ 
ity frequency pull ( 5oj cs = g 2 ^sx m) and photon lasing 
(Sk s = —2g 2 $tXm) [4, 74, 75]. At the phonon reso¬ 
nance (A m = 0), the Stokes linewidth equals k s + Sk s = 
(1 — C ) k s with C the same temporal cooperativity as 
above. Therefore, the threshold for Brillouin lasing is 
also C si. First realized in fibers [76], this case was re¬ 
cently also studied in CaF 2 resonators [14], silica disks [8] 
and chalcogenide waveguides [15]. Such lasers are known 
for their excellent spectral purity [75, 77] and received 
attention for quantum-limited amplification [16]. Since 
(n m ) = C(n s ) <C (n s ), there are far fewer phonons 
than Stokes photons in the cavity in this situation. The 


system enters the strong coupling regime as soon as 

Further, if the photon-phonon coupling rate is suffi¬ 
ciently strong, (27) simplifies to b = g 2 b. An identical 

derivation yields b = — g 2 b if the Stokes wave is consid¬ 
ered undepleted. Therefore, a red-detuned probe pro¬ 
duces entangled photon-phonon pair generation (b(t) oc 
e ±9t ), whereas a blue-detuned probe produces Rabi flop¬ 
ping between photons and phonons (b(t) oc e ±lgt ) [2]. A 
situation of equally strong optical and mechanical damp¬ 
ing (k s « K m ) invalidates the above weak coupling treat¬ 
ment even for small g. However, this is not sufficient to 
see strong coupling behavior. From the general solution 
(see Appendices), this requires g » k s + K m . Indeed, in 
the strong coupling regime the hybridized photon-phonon 
polariton sees half the optical and half the mechanical 
damping [12]. Therefore, the state-swap frequency - 
must be high compared to the average decay rate 
to observe an actual Rabi swap before the population 
decreases by 1 /e. 

By comparing (26) and ( 8 ), we prove an analogy be¬ 
tween spatial and temporal optomechanical effects (fig.3): 
the above cavity-based discussion still largely holds for 
guided-wave optomechanics with the mapping go ^ go, 
C h4 C, k s h 4 a s , K m h 4 a m and n <f>. In case of all 
co-propagating waves, and in the absence of cascading 
[30, 63] and noise [46, 47], the mapping of cavity optome¬ 
chanics onto a Brillouin-active waveguide in steady-state 
is an exact equivalence. However, when for instance one 
of the particles counter-propagates, such as the Stokes 
photons in backward scattering, important differences 
arise that have no equivalent in cavity optomechanics. 
Indeed, as proven in section III, information regarding 
the propagation direction of the waves disappears in the 
mean-field transition. Instead comparing (27) and (29), 
much can still be learned by instead mapping < 7 g H> — g 2 
and k s h 4 —a s . Note that this particular difference dis¬ 
appears if the counter-propagating particle species is un¬ 
depleted: then it vanishes from the dynamics and the 
situation is identical to the co-propagating case. 

Thus, guided-wave weak coupling requires g <C 
lets =F a m\ with g = go^/Wp the spatial coupling rate (see 
Appendices). Under weak coupling, there are two cases 
depending on the relative photon and phonon propaga¬ 
tion losses. We have touched upon these subcases at the 
end of section II and briefly consider them again here to 
show the similarity with cavity-optomechanical effects. 
First, when the phonons propagate far (cc m a s ), the 
photonic loss a s barely changes. However, the phononic 
response can then drastically change to y ” 1 — g 2 xt, which 
includes a shift in both the phononic propagation loss 
(<5a m = — 2 g 2 ?Ry*) and group velocity (oc SrXg); i.e. 
traveling-phonon amplification and light-induced slowing 
down of sound. In section V, we show that this unob¬ 
served regime can be achieved in existing systems. 

Second, when the Stokes photons propagate far (a s 
a m ), the phononic loss a m barely changes. However, 
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the photonic response can then drastically change to 
X^ 1 — g 2 Xm- Hence, we are back in the conventional 
domain of phonon-assisted amplification of traveling pho¬ 
tons (<5a s = — 2g 2 $tXm) and sound-induced slowing down 
of light (oc [38]- At resonance (A m = 0), the Stokes 
propagation loss is (1 — C)a s as in (9). 

If the coupling is sufficiently strong compared to the 
propagation losses (g a s +a m ), (29) simplifies to d 2 b = 
±g 2 b (see Appendices). In the forward (+) case, and with 
boundary condition 6(0) = 0, this yields 


phonon losses, so effectively lower C, in the cavity than 
in the waveguide. This directly shows that 

C > C (34) 

given (33), k 7 > t> 7 a 7 and C = . Clearly, the guided- 

wave cooperativity exceeds the cavity-based cooperativ- 
ity since the cavity has additional dissipation (e.g. cou¬ 
pling and bending losses). Finally, we define a gain coef¬ 
ficient Q for a cavity in analogy to (11) 


a s (z) = a s (0) cosh gz (30) 

b(z) = — io|(0) sinlipi; 


such that 4> s (z) = <f> s (0) cosh 2 gz and <f> rn (;z) = 

$ s (0) sinh 2 gz. Therefore, <& s (z) — $ m (2:) = <f> s (0) and 
<9 2 <f> s = d z <& m as required by the Manley-Rowe relations 
(4). In the backward (—) case, with L the waveguide 
length and boundary condition 6(0) = 0, the evolution 
along the waveguide has no cavity-optomechanics coun¬ 
terpart. Specifically, we retrieve 


a s (z) 

b(z) 


-—rcos gz 

cos gL 

.a\{L) . „ 

—i -— sm gz 

cos gL 


(31) 


such that <& s (z) = C o° 2 ( cos 2 gz and $ m (z) = 

cos^x sin2 9 Z - Therefore, $ s (z) + $ m (z) = and 

—(?,<!> s = 9 2 ( I> m as required by Manley-Rowe (4). The 
system has an instability at gL = 1] ■ which is reached 
before a full state swap between light and sound can be 
completed. This situation is called “contraflow Hermi- 
tian coupling” in classifications of coupled-mode interac¬ 
tions [78, 79]. In case of anti-Stokes (instead of Stokes) 
seeding in the strong coupling regime, an identical deriva¬ 
tion leads to d 2 b = —g 2 b - which produces the same 
Rabi oscillations for forward and backward scattering. 
Although familiar in resonators [2], these strong-coupling 
effects have not yet been observed in the field of guided- 
wave Brillouin scattering; see section V for the prospects. 

We conclude this section by analyzing the relation be¬ 
tween the guided-wave and cavity-based cooperativities 
(C and C ) and by introducing a gain coefficient (Q) for 
an optomechanical cavity. Note that the temporal coop¬ 
erativity 

4o 2 

C = —^ (32) 


is the ratio between the roundtrip gain and loss: inserting 
g 2 = 5o n p’ n p = and (21) in (32) indeed leads to 


C = 


Q Ly, V/n 0: n 


'"s ur 


QP V L v n 


k s T s 


(33) 


with P p the intracavity pump power and a natu¬ 

rally appearing correction factor that allows for higher 


4g 0 2 

hu]pK m 


(35) 


which characterizes the temporal exponential build-up of 
the Stokes when the phonons are heavily damped. The 
gain coefficients therefore obey 

Q > —Q (36) 

given n m > and (21). Hence, the guided-wave 

and cavity-based optomechanical figures of merit are now 
conceptually similar and the relations between each of 
them were given in (1), (21), (34) and (36). 


V. PROSPECTS 

In this section, we first give a couple of examples of how 
the Q O go connection (1) can be implemented - includ¬ 
ing several systems in which it can be tested empirically. 
Next, we move on to the prospects for observing new 
regimes of guided-wave optomechanics, simultaneously il¬ 
lustrating the application of our framework. Finally, we 
briefly discuss potential extensions to the minimal model 
(2) of traveling-wave Brillouin scattering. 

Table I presents four implementations of the conver¬ 
sion from the gain coefficient Q to the vacuum coupling 
rate go {Q —s- go) and four in reverse (Q go)- The sys¬ 
tems range from silicon nanowires and dual-web fibers 
to ultracold atom clouds and GaAs disks. In five cases, 
such as for silicon nanowires, the conversion can clearly 
be tested empirically by measuring Q and go through in¬ 
dependent, established methods [2, 23]. In three cases, 
the conversion is hypothetical but still allows for com¬ 
parison of the photon-phonon interaction strengths. For 
instance, an ultracold atom cloud in a Fabry-Perot cav¬ 
ity [55] has no obvious traveling-wave equivalent. Never¬ 
theless, its hypothetical waveguide partner would have a 
large gain coefficient of ~ 10 8 W -1 m -1 - which compares 
favorably to optomechanical waveguides realized to date. 

So far, all Brillouin-active waveguides had far lower 
phonon than photon propagation lengths (a" 1 a^ 1 ). 

Cavity-optomechanical systems, by contrast, more often 
than not had far lower phonon than photon damping 
rates (n m k s ) [2]. Only uniquely high-optical-quality 
[8, 14, 15, 74, 83, 84] systems succeed at reversing the 
latter hierarchy (/c m k s ). The common reversal of this 
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Q [Whn- 1 ] 

W 

£[H*] 

tHHz] 

Q m [ ] 

L [pm] 

«g [“I 

A [pm] 

Silicon nanowire [23, 24] 

10 4 


1.5-10 6 

V L [Mm] 

10 10 

10 3 

- 

4.6 

1.55 

Silica standard fiber [4] 

1 


70 

\/ L [cm] 

10 10 

500 

- 

1.45 

1.55 

Silica dual-web fiber [21] 

4- 10 6 


3-10 3 
\/ L [cm] 

6 • 10 6 

4- 10 4 

- 

1.7 

1.55 

Chalcogenide rib [80, 81] 

3- 10 2 


7-10 5 

V L H 

8 • 10 9 

230 

- 

2.6 

1.55 

Silica microtoroid [12] 

600 

<— 

3 • 10 3 

8 ■ 10 7 

2 • 10 4 

97 

1.45 

0.78 

Silicon optomechanical crystal [73] 

4- 10 4 


6 • 10 5 

6 • 10 9 

2 • 10 3 

5 

5 

1.55 

Rb ultracold atom cloud [55] 

10 s 


6 • 10 5 

4- 10 4 

42 

400 

1 

0.78 

GaAs optomechanical disk [82] 

5- 10 4 

<- 

3 • 10 5 

10 9 

2 • 10 3 

8 

4 

1.55 


TABLE I. Translation between waveguides and cavities. The table contains four conversions from a gain coefficient 
to a vacuum coupling rate and four in reverse. The right five columns contain the parameters necessary for the conversion 
through formula (1). These are order-of-magnitude estimates. In some cases, indicated by a *, the conversion can be empirically 
tested. In other situations, the conversion is hypothetical: e.g. an ultracold atom cloud in a Fabry-Perot cavity has no obvious 
guided-wave equivalent. 


damping hierarchy (going from waveguides to cavities) 
stems from the small phonon group velocities (v m <C v a ). 

The question naturally arises if waveguides with larger 
phonon than photon propagation length (a" 1 a" 1 ) 

can be made, while still keeping high cooperativities 
C = ~ 1. Currently, the largest phonon decay 

lengths are of the order a" 1 ~ 100 ;im in backward Bril- 
louin scattering [23, 65]. To realize larger phonon propa¬ 
gation lengths, one must look for waveguides with large 
acoustic group velocities v m and small linewidths K m . 
Thus, one promising approach uses low-frequency flex¬ 
ural modes (Q oc q 2 ) in backward mode (large q ) at low 
temperatures (large Q m ). Indeed, then we have both 
large v m oc q and small K m = ^ ~ Hz = 1 kHz 
[12, 85, 86]. Since a" 1 = ^ in a waveguide (where 
there is only propagation loss), we find that decay lengths 
up to a" 1 ~ 10 m are feasible given v m ~ 10 4 m/s and 
K m ~ 1 kHz. Such a small phonon propagation loss would 
strongly boost the cooperativity C, which could compen¬ 
sate for a potentially lower coupling rate g in backward 
mode. Clearly, nothing intrinsically forbids amplification 
of traveling phonons in systems such as the dual-web fiber 
[21] - where a' 1 ~ 10 cm. Besides its scientific inter¬ 
est, such a traveling-phonon amplifier may be useful in 
phonon networks [13, 48, 49, 87, 88]. 

Next, we look into achieving the strong coupling regime 
in the typical situation of high acoustic loss (a m a s ). 
To see traveling-wave Rabi flopping, entangled photon- 
phonon pair production or contraflow Hermitian coupling 
(see section IV), one must obtain g > a m or equivalently 
C > yp 1 . In optical fibers, in backward mode and given 

a” 1 ~ 10km and a” 1 ~ 100pm, this requires C > 10 8 . 
This necessitates an unrealistic continuous-wave pump 
power of P p > 10 8 ^* = lOkW with Q ~ 1 W _1 m _1 . 


In contrast, silicon chips can produce significantly 
lower ^ ratios and therefore ease the condition on C 

OL s 

for strong coupling. One can expect phonon propaga¬ 
tion distances up to a” 1 ~ 1 mm, as these are readily 
achieved in surface-acoustic-wave devices [89]. Together 
with a" 1 ~ 1cm [23], this yields C > 10 as the strong 
coupling condition, which requires a reasonable pump 
power of P p > 10 = 100 mW with Q ~ 10 4 W _1 m _1 

[23, 24]. Indeed, current nanoscale silicon systems have 
already demonstrated C « 2 [24, 63]. Hence, taking 
into account the rapid progress in state-of-the-art devices 
[1, 23, 24, 63], we expect demonstrations of traveling- 
phonon amplification and guided-wave strong coupling 
in the coming years. Such observations would open up 
entirely new realms of optomechanics. 

Finally, this work can be extended on several fronts. 
First, the mean-field transition can be applied to the 
noise models of [46, 47]. Second, the regime of “nonlin¬ 
ear quantum optomechanics” [2, 90-92] should be trans¬ 
ferred to waveguides. This requires that strong coupling 
is reached for merely one pump photon (<f> p i—> Is -1 ): 
9 = 9o\/% = go > « s + a m . As a m > a s usu¬ 
ally, traveling-wave nonlinear quantum optomechanics is 
achieved when go > a m . Third, the coupling between the 
phononic mode and the thermal bath [46, 47] must be 
treated carefully to obtain truly quantum-coherent [12] 
coupling. And fourth, we focused mainly on the dynam¬ 
ics that optomechanical waveguides and cavities have in 
common, but wisdom may be found in the differences as 
well. We gave the example of contraflow strong coupling 
in section IV. In addition, the cavity has input fluxes 
that have no equivalent in a typical guided-wave set-up, 
while the waveguide can display spatiotemporal effects 
(both d z and dt in (2)) that are absent in a cavity (only 
dt in (20)). On top of this, the cavity breaks the symme- 
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try between Stokes and anti-Stokes scattering, whereas 
this symmetry prevents exponential build-up of noise in 
low-dispersion forward intra-modal scattering [47]. It has 
also yet to be determined whether different cavity dy¬ 
namics results in the medium-finesse case, as section III 
was limited to a low or high finesse. 

With slight modifications, (2) also captures Raman ef¬ 
fects [4, 6, 93]. For instance, the phonon frequency is 
much larger so an optical phase-mismatch can arise. Fur¬ 
ther, the thousandfold higher optical phonon frequency 
puts most Raman modes in the ground state at room 
temperature. This breaks the symmetry between Stokes 
and anti-Stokes scattering, such that thermally seeded 
exponential build-up may be seen even in the forward 
intra-modal case [47]. Some parameters introduced here 
- such as the spatial vacuum coupling go - lose elegance 
in the Raman case: the impossibility of significant opti¬ 
cal phonon transport undermines their symmetric defi¬ 
nitions. Mathematically, expressions such as (5) diverge 
as v m and a” 1 11 vanish. Then one must resort to the 
broken-symmetry Raman gain coefficient, which never¬ 
theless obeys (1). Therefore, the core of this work also 
applies to guided-wave [4, 94-96], cavity-based [97-100] 
and surface-enhanced Raman scattering [57]. 

VI. CONCLUSION 

In conclusion, we unveiled a connection between 
Brillouin-active waveguides and optomechanical cavities. 


The link between the Brillouin gain coefficient Q and the 
zero-point coupling rate go was derived in a platform- 
independent way. As illustrated for silicon nanowires 
and ultracold atom clouds, it significantly expands the 
variety of systems whose photon-phonon interaction effi¬ 
ciency can be compared. Through the mean-field tran¬ 
sition, we derived the dynamics of optomechanical cavi¬ 
ties from that of Brillouin-active waveguides. We framed 
the behavior of both systems in terms of cooperativi- 
ties and vacuum coupling rates. Next, we showed that 
phenomena familiar from cavity optomechanics all have 
guided-wave partners, but not the other way around. 
The broader theory predicts that several novel regimes, 
such as guided-wave strong coupling, will be accessible 
in state-of-the-art systems in the coming years. Hence, 
we showed that Brillouin scattering and cavity optome¬ 
chanics are subsets of a larger realm of photon-phonon 
interaction. 
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Appendix A: Link to cavity Hamiltonian 


Appendix B: Manley-Rowe relations 


With the mean-field transition derived in the main 
text, we take a step beyond the Q go link. As we show 
in this section, the mean-field equations are equivalent to 
the cavity Langevin equations in the resolved-sideband 
limit (k -C fi m ). In the case of coupling between one 
mechanical and one optical resonance (fig.2a), the usual 
theory [2] starts from the Hamiltonian 

H = huj c a)a + + Hint 


with 


Hint = hgoata (b + 6 f ) 


the interaction Hamiltonian, x = xzpf (b A U j the me¬ 
chanical oscillator’s position, Xzpf the zero-point motion, 
a and b ladder operators for the optical and mechanical 
oscillator and go = xzpf the zero-point coupling rate. 
When the pump is undepleted, the interaction Hamil¬ 
tonian can be linearized: a = a + 5a with 5a a small 
fluctuation. Then we have 

Hin t n ) = hg 0 a (5a + 5at) (b + 


Using the equation of motion a = — 4[a, H] and the 
commutator [o, fit] = 1 (the same for b), this linearized 
Hamiltonian leads straight to the coupled equations [2] 

5a = — — iAj 5a — igoa (b + b ^ 

b = — A — iCl n ^j b — igoa (5a + 5a?) 

with A = oj p — lo c . Next, we consider a blue-detuned 
pump in the resolved-sideband regime (k -C fi m ). Then 
we can write the ladder operators as 5a —> a s e int and 
b —> be~ int , with o s and b now slowly-varying. We neglect 
the 6-term in the optical equation and the <5a-term in the 
mechanical equation because they are off-resonant. This 
is the rotating-wave approximation, which corresponds 
to the classical slowly-varying envelope approximation [3, 
4]. Hence, the above equations reduce to 


a s = -igoutf - x s (Al) 

b = -ig 0 aal - x^b 


and we find that equations (Al) are identical to equations 
(16) given a s H* a s and b K> b. Remarkably, the equiva¬ 
lence holds even though the pump and Stokes could be 
counter-propagating or in different optical modes. In the 
unresolved-sideband limit (f2 m <C k), anti-Stokes gener¬ 
ation and cascading must be added for forward intra- 
modal, but not necessarily for backward or inter-modal 
Brillouin scattering. Indeed, comb generation is usually 
not accessible by backward or inter-modal coupling be¬ 
cause of the phase-mismatch (fig. 2). This assumption 
can be violated in Fabry-Perot cavities [101] or when the 
first-order Stokes becomes sufficiently strong to pump a 
second-order Stokes wave [9]. 


In this section, we prove that the Manley-Rowe rela¬ 
tions guarantee the existence of a single real, positive 
photon-phonon coupling coefficient in waveguides (go) 
and in cavities (go)- In waveguides, the Manley-Rowe re¬ 
lations are formulated at the level of photon and phonon 
fluxes $. In cavities, they are written down in terms of 
the total photon and phonon numbers n. 


a. Manley-Rowe in waveguides 


A Brillouin-active waveguide in steady-state (dt —> 0) 
obeys (see (2)) 

H up — iK mop a s b x p ftp 
±d z a s = -ikmostfap - x“ 1 a s (Bl) 

d z b = —iR ora a\a p - x^b 

with arbitrary normalizations of the pump, Stokes and 
acoustic envelope such that generally k mop i=- k mos i=- 
kom are different complex numbers. Using d z (ft^ft) = 
(d z a))a + a^dza), we find 

<9,<f>p = -Qp-Pp + (iR* mop a\b ] a v + h.c.) 

±<9 z $s = -a s <P s - (*Kmos4 &ta P + h.c.) (B2) 

rn — A h.C.) 

Suppose now that the envelopes are flux-normalized such 
that <f> p = ftfft p , <f> s = a|ft s and <!> tr , = b^b correspond 
to the number of pump photons, Stokes photons and 
phonons passing through a cross-section of the waveg¬ 
uide per second. Then we demand that, in the lossless 
case (cq = 0), the rate of pump photon destruction equals 
the rate of Stokes photon and phonon creation 


~d z % = ±ch$ s = (B3) 

These are the Manley-Rowe relations [3, 79] for a Bril¬ 
louin waveguide. We deduce from (B2) and (B3) that 

^mop — Anns — Korn (B4) 

This proves the existence of a single coupling coefficient 
that captures all reversible optical forces and scattering. 
Note that (B4) also guarantees power-conservation since 

d z (HuipQp ± hoj s ^ s + hfl$ ul ) = 0 

leads with (B2) in the lossless case to 


K*'p K|nop A ^s^mos A — 0 (B5) 

which is true given (B4) and co p = u; s + H. Next, we show 
that this coefficient (B4) can be taken real and positive 
without loss of generality. Renormalizing the envelopes 
to Cpftp, c s a s and c m b yields new coupling coefficients 


Cs Cm 


C, 
CpC s 


—lA-ftom (B6) 


CpC m 
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as can be seen from (Bl). Suppose that K om = goc lip is 
complex with go real and positive. Then we take c p = 
c a = c m = e _ *A Using (B4) and (B6), it follows that the 
renormalized coupling coefficients are real and positive: 

i^mop — it rn os — itom — 90 (B7) 

This unique coupling coefficient quantifies the coupling 
strength between a single photon and a single phonon 
propagating along a waveguide. Indeed, suppose that 
a p = a s = i) 1 s” 1 / 2 such that <f> p = $ s = <f> m 1-4 1 s -1 
at a certain point along the waveguide. In the lossless 
case, (B2) then becomes 

d 2 $ p = -250 

±d z <$> 3 = 2g 0 (B8) 

d z *& m = 2^0 


So 2go gives the rate (per meter) at which the pump flux 
decreases and the Stokes and phonon flux increase at a 
point along waveguide through which one pump photon, 
one Stokes photon and one phonon are passing. 

The waveguide coupling coefficient go can also be in¬ 
terpreted in terms of a zero-point motion. As shown in 
(14), we have 


9o 


L 


v p v a v m 


-go 


(B9) 


For forward intra-modal scattering (v p = v s = v g ) 


go = xzpf 


dujp 

dx 


kp 


(BIO) 


is defined in terms of the zero-point motion and the cavity 
frequency pull at fixed wavevector [2]. Combining (B9), 
(BIO) and (D9), we obtain 


9 0 = ~ 


-£ZPF 

C 


dn eS 


dx 


= ~ %ZPF 


dk p 
dx 


(Bll) 


with 


*zpf = = \j pppp < B12 > 

the waveguide “zero-point motion” and m e ff the effective 
mass per unit length. Indeed, a waveguide section of 
length L contains on average (n m ) = ^p-($ m ) phonons 
with ($ m ) the mean phonon flux. As fluxes - instead of 
numbers - are the fundamental quantities in waveguides, 
the zero-point motion is corrected by precisely a factor 

\fE in ( Bl2 )- 

Often the optical envelopes are power-normalized and 
the acoustic envelope displacement-normalized. Starting 
from flux-normalized envelopes, one can switch to such 
normalizations through 


Cn — 


Cs = 


Cm — 


2 hfi m 

^AlfAri 


= 2izPF 


(B13) 

with k e s the effective stiffness per unit length and by 
applying (B6). 


b. Manley-Rowe in cavities 


Here, we apply the discussion of the previous section to 
the mean-field cavity equations. With arbitrary envelope 
normalizations and without input, equations (13) are 

Up — in mop a s b y p u p 

a s — in rnosb u p y s u s (B14) 

b — iK, om d a cip X m b 

with generally K mop ^ K mos J=- K om . Applying ^ (a f u) = 
(g^u^)u + a^(gju) to (B14), we find 

^ n p = —KpUp A (iK* wp alPa p A h.c.) 

^ n s = —K s n s - (iK mos alb f a p A h.c.) (B15) 

— n m n m (^in om cpb u p A h.c.^ 

Suppose now that the envelopes are number-normalized 
such that n p = aj,u p , n s = aju s and n m = Pb corre¬ 
spond to the number of pump photons, Stokes photons 
and phonons in the cavity. We demand that, in the loss¬ 
less case (Kj = 0), the rate of pump photon destruction 
equals the rate of Stokes photon and phonon creation 


-iip =h s = h m (B16) 

These are the Manley-Rowe equations for an optome¬ 
chanical cavity. We deduce from (B15) and (B16) that 

^mop = ^mos — ^om (B17) 

This proves the existence of a single coupling coefficient 
that captures all conservative optical forces and scatter¬ 
ing. Note that (B17) also guarantees energy-conservation 
since 

— (HaipUp A hco s n a A = 0 

leads with (B15) in the lossless case to 

W p K mop A LO a K rnoa A fl/Com — 0 (BIS) 

which holds given (B17) and uj p = oj s A fi. As in the pre¬ 
vious section, one can show that this coupling coefficient 
can be chosen real and positive. This unique coupling 
coefficient must then be the well-known go- It quantifies 
the interaction strength between a single photon and a 
single phonon trapped in a cavity. Indeed, suppose that 
a P = a s = b 1-4 1 such that n p = n a = n m 4 1 at a 
certain point in time. In the lossless case, (B15) then 
becomes 


11 p = -2ffo 

n s = 2<7o (B19) 

n m = 2g 0 
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So 2go gives the rate (per second) at which the number 
of pump photons decreases and the number of Stokes 
photons and phonons increases when there is one pump 
photon, one Stokes photon and one phonon in the cavity. 

Often the optical envelopes are energy-normalized and 
the acoustic envelope displacement-normalized. Start¬ 
ing from number-normalized envelopes, one can switch 
to such normalizations through 


= s/hui p c s = \JttiJs 



2xzpf 


(B20) 

with xzpf the zero-point motion and by applying (B6). 


Appendix C: Mean-field approximation 

c. Justification of fg = fg 

We denote f(z , t) and g(z , t) two operators that vary 
slowly on a lengthscale L. The mean-field operators 
are defined as f(t) = f(z,t)dz. Clearly, when 

f(z,t) = f(0,t) and g(z,t) = g(0,t) are constants then 
fg(t) = f(0,t)g(0,t) = f(t)g(t). Let us assume now 
that the amplitudes vary slowly enough such that they 
can be Taylor-expanded as f(z,t) = f(0,t) + f'z with 
f = <9-/(0, t ) and the same for g. Then we have 


d. From independent full-vectorial definitions 

Here, we derive equation (Dl) from the full-vectorial 
definitions of Q and go - specializing to intra-modal for¬ 
ward scattering. We focus on the moving boundary con¬ 
tribution. From the perturbation theory of Maxwell’s 
equations with respect to moving boundaries [102], the 
cavity frequency shift can be expressed as 

doj c _ L 0 P § dA (u • n) (Ae|E|| | 2 - Ae -1 |D ± | 2 ) 
dx ~ 2 JdVe |E| 2 

with u the normalized (max(|u|) = 1) acoustic field, n 
the unit normal pointing from material 1 to material 2, 
Ae = ei — £2 and Ae -1 = e/ 1 — e/ 1 . The upper integral is 
over the entire surface area of the cavity, the lower inte¬ 
gral across the cavity volume. Further, E|| is the electric 
field parallel to the boundary and the displacement 
field perpendicular to the boundary. For a longitudinally 
invariant cavity, the surface integral can be reduced to a 
line integral and the volume integral to a surface integral: 

du c Up §dl{ u.n)(Ae|E || i 2 -A e - 1 |D ± | 2 ) 
dx 2 fdAe |E| 2 1 ’ 

Further, the gain coefficient Q is given by [23, 25, 103] 

^ =Wp Jh l(f ’ u)|2 (D3) 


(/(O)/- I /'<) 

<1 ) |#)i + fi'(0)yj 

where we dropped the time-dependence. Thus, we have 

Jg = f(0)g(0) + (f'g(O) + f(O)g') \ + fg’^ 
Similarly, 

Jg = \ £ mm + if 9(0) + /(O)fl')« + fg’z 2 ) d^ 
= f(0)g(0) + (fg( 0) + mg’) § + /Vy 


with f the power-normalized optical force density and 
(f, u ) = / f • u dA. Note that fe e fi is the effective stiffness 
per unit length. In the case of radiation pressure forces 
f rp we have [103] 

f rp = ^ (Ae|e|| | 2 - Ae _1 |dj_| 2 ) n5(r - r bo undary) 

with <5(r — rboundary) a spatial delta-distribution at the 
waveguide boundaries. The fields e and d are power- 
normalized. Here we already assumed that the Stokes 
and pump field profiles are nearly identical, which holds 
for intra-modal SBS given the small frequency shifts. 
Hence, we get 

(f rp ,u) = | j) dl (u • n) (Ae|e||| 2 - Ae _1 |d_L| 2 ) (D4) 

Additionally, the guided optical power P is given by 


Therefore fg — fg = f'g' y/ which can be neglected if L 
is sufficiently small compared to the lengthscale on which 
f(z , t) and g(z, t) vary. 


Appendix D: Alternative proofs of the Q -O- go link 

In this section, we describe two other approaches to 
derive the link 


2 2 v^p) 

%=v z^iT- 



(Dl) 


P = — (E, eE) = 

2 ' ' 2 


eL4e|E| 2 


(D5) 


Combining equations (D2), (D4) and (D5), we find 
d w c _ v s uj p 
dx ~ 2 { vp ' ' 

A similar derivation can be done for the strained bulk, 
so we have 

duj c VgLOp 


dx 

( f - u ) = 


<f» 


2 

2 du) c 
v g ujp dx 


(D6) 
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with f = f rp + f es and f es the electrostrictive force density. 
Substituting equation (D6) in (D3) yields 


Substitution of (DIO) in (Dll) with xzpf = \J 2m 
results in 


2 Qm 

UlpVjkeff V dx ) 


(D7) 


4TQ m 2 

fiWpVjQm 9 ° 


Finally, we use the definition of the zero-point coupling 
rate go = aJzPF^f and the zero-point motion xzpf = 


2m e ffLO n 


with TO e ff the effective mass per unit length. 


Inserting these in (D7) yields 


Q 


2 Qn 


WpP|fc e ff 
4 L 


2m e ffFf2 m 2 

-7- 9o 


Q 


9o 


1 {huip) fi m v\ 


(D8) 


or the other way around 


2 2 y^v>) 

9 « = v *^r~ 



(D12) 


This proof only holds for forward intra-modal scattering 
- whereas the mean-field transition applies to backward 
and inter-modal scattering as well. 


Appendix E: Full solution of guided-wave evolution 


and (D8) is identical to (Dl). In this derivation, we 
started from full-vectorial definitions that are only valid 
for intra-modal forward scattering. In contrast, the 
mean-field transition shows that this result remains true 
with v g —> yfv p v s for inter-modal coupling. 

e. From independent derivative definitions 


In this section, we give the full solution of the traveling- 
wave spatial dynamics (8) in the constant pump ap¬ 
proximation ($p(2;) = 4>p(0)) where the pump is strong 
enough to be treated classically. From this solution, one 
can derive the regimes treated in section IV as limiting 
cases. In addition, this solution can directly be mapped 
on the cavity-based temporal dynamics (26). Thus, we 
start from 


The cavity resonance condition is k p L = 2nm with m 
an integer. Given k p = ^ 1 ’" aff and c the speed of light, 
this implies that 


d z ci s — T Xs (El) 

d z h = -ig 0 ala p - x^ h 


diOp 

dx . n e ff dx , 

<vp u-p 

This can be recast in terms of the index sensitivity at 
fixed frequency by 


which immediately yields 

(d z + Xm*) {dz ± Xs _1 ) a s(z) = ±9 2 a s {z) (E2) 

where ± stands for forward (+) and backward (—) scat¬ 
tering. Inserting the ansatz a s (z) oc e JZ leads to 


dn e ff 
dx fep 


Heff dn e ff 
n„ dx 


with u p h = the phase velocity and n g = the group 
index. Thus we have 


duj p 

dx 


hp 


uj p dn e ff 


fip _ 

n„ dx 


(D9) 


The cavity frequency pull must be calculated at fixed 

wavevector (go = £zpf yyy ), so this yields 

I fc D 


dn e ff 


dx 


2 ( W P 

= 9o ■'Kzpf — 


(DIO) 


7 2 + (Xm* ± Xs *) 7 ± (x m *x s 1 - 9 2 ) = 0 (E3) 

Its solution is 

71 = \ | - (Xm* ± Xs _1 ) + \J (Xm*TXs -1 )' +45 2 | 

(E4) 

72 = \ | - (Xm* ± Xs -1 ) - \/ (Xm*TX s _1 )“ +4^ 2 | 

Given these two roots, one can determine the exact evo¬ 
lution of the photon-phonon fields along the waveguide if 
the correct boundary conditions are known. The bound¬ 
ary condition 6(0) = 0 and fixed probe flux $5(0) = 
a|(0)o s (0) are appropriate for forward scattering such 
that 


Previously [23], we showed that 


Q = 2cj, 


Qm [1 dTl e ff 


&eff \ C dx 


(Dll) 


a s (z) = Qs ^ {(72 + X s *) e7lZ - ( 7 i + X s *) e722 } 


72 - 7 i 


(E5) 


b\z) = i a 4° ) (72 +^ 1 )( 7l+ ^ 1 ) { e»‘ - e»’} 


g 


72 - 71 
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The backward case (fixed <I> S (L) with L the waveguide 
length) can be solved similarly. This full solution con¬ 
tains the important regimes discussed in section IV. For 
instance, in the strong coupling regime (g a s + a m ) 
and at resonance (Aj = 0) one can show that 


71 

72 


(*m + Q; s large g _ 

- —+9 —> 9 


~9 


-g 


(E6) 


actual photon-phonon pair generation can be seen. The 
photons and phonons indeed equally share the total prop¬ 
agation loss a m + a s in this regime, as in cavity settings 
[12]. The spatial evolution (E5) then becomes identical 
to (30). The weak coupling regimes of stimulated photon 
(cti s <C a m ) and phonon (a m <C a s ) emission can equally 
well be obtained from the full solution (E5). This solu¬ 
tion also contains acoustic (a s <C a m ) [65] and optical 
{a m <C a s ) build-up effects. 


Therefore, the spatial coupling rate g must overcome the 
average photonic and phononic propagation loss before 



